How the quark self-energy affects the color-superconducting gap 
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We consider color superconductivity with two flavors of massless quarks which form Cooper pairs 
with total spin zero. We solve the gap equation for the color-superconducting gap parameter to 
subleading order in the QCD coupling constant g at zero temperature. At this order in g, there is also 
- - , a previously neglected contribution from the real part of the quark self-energy to the gap equation. 

' _ _| ' Including this contribution leads to a reduction of the color-superconducting gap parameter (/jq by a 

^ \ factor b'o = exp [ - (tt^ -I- 4)/8] ~ 0.177. On the other hand, the BCS relation Tc ~ 0.57(^o between 

■ and the transition temperature Tc is shown to remain valid after taking into account corrections 

from the quark self-energy. The resulting value for Tc confirms a result obtained previously with a 
\ different method. 

O: 

! I. INTRODUCTION 

^ ' Quantum chromodynamics (QCD) is the fundamental theory of the strong interaction. In strongly interacting mat- 
\^ ^ ter at large density or, equivalently, large quark chemical potential /i, asymptotic freedom |Q implies that single-gluon 
T-H . exchange becomes the dominant interaction between quarks. Single-gluon exchange is attractive in the color- antitriplet 
' channel By Cooper's theorem [Q, any attractive interaction destabilizes the Fermi surface and, at sufficiently 
, small temperature T, leads to the condensation of Cooper pairs. If the Cooper pair condensate carries charge quantum 
\ numbers of a local gauge symmetry, the Meissner effect leads to superconductivity. Strongly interacting matter, where 
quark Cooper pairs carry color charge, becomes a color superconductor. In a superconductor, exciting particle-hole 
pairs costs at least an energy amount 2(j)Q, where (j)^ is the value of the superconducting gap parameter at the Fermi 
surface for T = 0. Its value can be computed from a gap equation derived in mean-field approximation, which, in 
I QCD, involves single-gluon exchange |^,|^. 

^ Schematically, this gap equation can be written in the form p] 



o 



K = .9' 



(1) 



For small values of the QCD coupling constant, g ^ 1, the solution is 

</>o = 2 6/iexp(^-^j [1 + 0(5)] • (2) 

The first term in Eq. (|^) contains two powers of the logarithm ln(/i/0o)- One logarithm is well-known from the gap 
equation in standard BCS theory where it arises from the integration over fermion momenta up to the Fermi 
surface. The other logarithm is special to theories with long-range interactions, like the exchange of almost static 
magnetic gluons in QCD Its origin is a collinear singularity when integrating over angles between quark and 

gluon momenta in the gap equation. The weak coupling solution (^ implies that the first term in brackets in Eq. (Q) 
is ~ It therefore dominates the right-hand side of Eq. (|l|). Together with the prefactor , it is of order 0(1) 

in the gap equation. The value of the coefficient Q determines the constant c in Eq. (||). 

The second term in Eq. (|l]) contains subleading contributions of order 0{g) to the gap equation, characterized by 
a single power of the logarithm ln(/i/0o) ^ ^/g- A part of these contributions arises from the exchange of non-static 
magnetic and static electric gluons [||. Both types of interactions are short-range: they are screened on a distance 
scale rn~^, where nig is the gluon mass; = N f g'^ ji^ / (Qir^) , Nf is the number of quark flavors. Consequently, the 
collinear logarithm characteristic for long-range interactions is absent, and one is left with the BCS logarithm. The 
coefficient P in Eq. (^ determines the constant b in Eq. (||). 

The third term in Eq. ([T]) summarizes sub-sublcading contributions of order 0{g^) with neither a collinear nor a BCS 
logarithm. It was argued in that at this order gauge-dependent terms enter the QCD gap equation. However, 
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the gap parameter is in principle an observable quantity, and thus gauge-independent. Therefore, one concludes that 
the mean-field approach cannot be used to compute sub-subleading contributions to the gap parameter. It was also 
shown that effects from the finite lifetime of quasiparticles in the Fermi sea influence the value of (fiQ at this 
order. In weak coupling, these contributions are suppressed by one power of g compared to the subleading terms and 
therefore constitute an order 0{g) correction to the prefactor 6, as indicated in Eq. (^. 
The value of the coefficient c was first computed by Son j^] , 



TT 



25 ' '~ 



(3) 



Son also gave an estimate for the constant b, 



with a constant bo of order 0(1), which could not be determined in the approach of Ref. 0. In 1^,^ the constant 
bo was computed by solving the QCD gap equation including non-static magnetic and static electric gluon exchange. 
The result is 

60 = 256^4 (^—j b'o, (5) 

with an undetermined constant b'o of order 0(1). In |^,||, where the only subleading contributions to the gap equation 
arise from non-static magnetic and static electric gluon exchange, 6q = 1. In principle, however, there could be other 
subleading contributions, which would change 6q to a value &q 7^ 1. 

At sufficiently large temperature, thermal random motion breaks up Cooper pairs and the superconducting con- 
densate melts. In Ref. ||^, it was shown that the temperature Tc for the transition between the normal and the 
superconducting phase is related to the zero-temperature gap at the Fermi surface in the same way as in BCS theory, 

T, ^ — ^0 - 0.57 (j)o , (6) 

TT 

where 7 ~ 0.577 is the Euler-Mascheroni constant. 

In Ref. ijl^ , Brown, Liu, and Ren calculated in a different approach with the result 



where 



c'l^expl-^^'l 0.177. (8) 



Furthermore, the authors of ||l^ assumed the validity of Eq. (^), and concluded that 

b'o = c'l , (9) 



as one readily checks with Eqs. (||)-(||), (Q), and (||). Physically, the difference between the approach of Refs. 
and that of [ p2| is that contributions from the quark self-energy were neglected in the former, but taken into account 
in the latter. If the above arguments are correct, one may therefore conclude that the quark self-energy constitutes a 
subleading correction to the gap equation and thus is responsible for the change of the value of b'o from 1 to c'l given 
by Eq. (||). The authors of also assert that there are no further subleading contributions that could alter the 
value of c'l . 

While manifestly gauge invariant, a disadvantage of the approach of ||l2| is that its range of applicability is restricted 
to the normal-conducting phase and thus can only determine the value of Tc, but not the value of the zero-temperature 
gap 00. Therefore, a relation between Tc and i^oi like Eq. (^, can in principle not be established within this approach. 
It is possible to derive such a relation with the help of the gap equation, as demonstrated in j5|, but contributions 
from the quark self-energy were neglected in obtaining the result Eq. (||). It is conceivable that Eq. (|^) changes, once 
these contributions are taken into account. Consequently, the validity of Eq. is not obvious. 

The aim of the present paper is twofold. On the one hand, we want to compute the contribution of the quark 
self-energy to the value of the constant b'o in the zero-temperature gap. On the other hand, we want to confirm the 
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result (|7|) for T^. To this end, it is necessary to first compute the value of the zero-temperature gap by directly solving 
the gap equation including the quark self-energy. Second, one has to prove that Eq. ^ remains valid in order to 
determine Tc, which then can be compared to the value obtained in [|2|. Our paper is organized as follows. In 
Section II we first clarify how the quark self-energy enters the gap equation. In Section III, the resulting gap equation 
is solved at zero temperature. In Section IV we determine Tc. We conclude in Section V with a summary of our 
results. 

Our convention for the metric tensor is g^'^ = diag{l, —1, —1, —1}. Our units are h — c = ks ~ 1- Four-vectors are 
denoted by capital letters, K = = (fco, k), and k = |k|, while k = k/fc. 



II. THE GAP EQUATION INCLUDING THE QUARK SELF-ENERGY 

In fermionic systems at non-zero density, it is advantageous to treat fermions and charge- conjugate fermions as 
independent degrees of freedom and to work in the so-called Nambu-Gorkov basis. In this basis, the full inverse 
fermion propagator is defined as lll| 



s- 







H ] 








^22 / 





S21 



'0-1 . ■ 



(10) 
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where 5**^11 is the propagator for free fermions, S''^22 the propagator for free charge- conjugate fermions. In momentum 
space and for /i ^ to. 



S\i{Q) = [i^Q, + /no)"' , ^°22(Q) = irQ, - /no)"' 



(11) 



where 7'^ are Dirac matrices. The four components of the fermion self-energy are denoted as , i,j — 1, 2. The 11 
component of the self-energy, Sn, is the standard one-loop self-energy for fermions; similarly, E22 is the self-energy for 
charge-conjugate fermions. The 21 component of the self-energy, S21, which was denoted in is the gap matrix 
in a superconductor, while S12 = 7oS2i7o- In the following, we somewhat imprecisely use the term "self-energy" only 
for the diagonal components En and £22- 

Inverting Eq. (nOf) one obtains the full fermion propagator S, with the diagonal components 



•S^ll + Sii — I]i2(5°22 



5*11 — 

describing the (normal) propagation of fermions, and 



S22) 'S21 



(12a) 



^2 



E22 — S2l(S'°ii -I- Sii) 



(12b) 



describing the (normal) propagation of charge-conjugate fermions. In superconductors, due to the presence of a 
fermion-fermion condensate one can always convert an incoming fermion into an outgoing charge-conjugate fermion 
and vice versa. Therefore, the full fermion propagator S also has off-diagonal components. 



S'12 — ^('5''^ii + Sii) "'"Si2S'22 , 
»S'21 = ^(>5''^22 + ^22) "'"S2lS'll , 



(12c) 
(12d) 



describing the anomalous propagation of fermions and charge-conjugate fermions. [Please note that our sign conven- 
tion for the self-energy differs from that in ||ll[], which leads to the difference between our Eqs. ( |l2[) and Eqs. (2.4) 
and (2.5) in Q.] 

Let us now consider a system of quarks interacting via one-gluon exchange. In mean-field approximation [T^ , the 
four components of the fermion self-energy in momentum space are computed as 



i,i = i,2 



(13) 



Here, is the gluon propagator, and is the diagonal Nambu-Gorkov matrix = diag(7'^TQ, —^^T^ ), Ta are the 
Gell-Mann matrices. We compute the self-energy in the imaginary-time formalism, i.e., T /V J^q = ^ / 'i/ i'^'^)^ j 
where n labels the fermionic Matsubara frequencies, w„ — {2n + l)TrT = iqo. For ij — 11, Eq. ( p^ becomes Eq. (2.7) 
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of 0, for ij = 21, we recover Eq. (2.6) of (however, due to our different sign convention, only up to an overall 
sign). 

A fully self-consistent treatment of the mean-field approximation requires to solve the coupled system of Eqs. ( [T^ ) 
and (13). The mean-field solution obtained in this way resums terms of infinite order in the coupling constant. 
However, because only a particular class of diagrams is taken into account (the so-called "rainbow" diagrams), 
such a solution is in general not gauge invariant. On the other hand, the quasi-particle properties encoded in the 
propagator, like their excitation spectrum, are physical observables and thus in principle gauge invariant. Indeed, 
a complete solution of the Schwinger-Dyson equations, as well as a perturbative expansion in powers of g, preserve 
gauge invariance. Nevertheless, as was discussed in the introduction, an expansion of the mean-field equation ( p^ ) 
for the color-superconducting gap matrix = E21 in powers of g is believed to be gauge invariant up to terms of 
subleading order [the first two terms in Eq. ([^)], and the gauge dependence only surfaces at sub-subleading order [the 
third term in Eq. (|l]), or the terms ~ 0(17) in Eq. (^)]. To preserve gauge invariance beyond subleading order, other 
diagrams than those of rainbow topology have to be added to Eq. (p^, or in other words, one has to go beyond the 
mean- field approximation to solve for the gap matrix £21- 

If one restricts the computation of the gap to subleading accuracy, however, the mean-field equation (|l3| ) should 
be sufficient to obtain a gau ge-invariant result. It is then mandatory to identify all terms that can contribute to 
subleading order . In Eq. ( [l3|) for S21, the term in square brackets becomes —^^^T^ S2i{Q)j'^Tb. With the exception 
of Ref. Ill, previous calculations of the QCD gap parameter neglected the terms En and E22 in 5*21, see Eq. 



(p^). A perturbative calculation of these self-energies, i.e., approximating f^S'((3)f^ 
analogously for ij = 22, and analytical continuation to real energies qo gives the result 



E«(g) = S°n(Q) = S022(Q) 



'lag [ qoln^- +i7r|go 
% 



l^TaS\i{Q)YTb, and 



(14) 




where = (37r/4)TO^. On the quasi-particle mass shell, go — ^q, and near the Fermi surface, ~ (j^o, the real part 
of the self-energy is of order g^c/iQ ln(/i/0o) ~ 3</'0) while the imaginary part is ^ g'^4'0 and thus down by a factor of g 
compared to the real part. 

In |l^, the real part of YP was neglected and the effect of the imaginary part on the magnitude of the color- 
superconducting gap was studied. It was found that a non- vanishing imaginary part leads to sub-subleading corrections 
[terms included in the third term '-^ a in Eq. (|l|)] and to corrections of order 0{g) to the prefactor of the gap, cf. Eq. 
Therefore, they are of the same order as terms that violate gauge invariance in the mean-field approximation 

Since the real part of the self-energy is parametrically larger than the imaginary part by one power of (/, we expect 
the former to contribute to subleading order, 0(g), to the gap equation, and therefore lead to a correction of order 
0(1) to the prefactor. As discussed in the introduction, this is precisely what the authors of [|l^ found, assuming the 
validity of Eq. (^. In the next section, we solve the gap equation including the quark self-energy and compute the 
value of 5q at zero temperature. In Section IV we then check the validity of Eq. (g). 

First note that, since the real part of the quark self-energy is expected to influence the value of 4>o only at subleading 
order in the gap equation, it is sufficient to approximate the value of En or E22 in the propagator in Eq. (|l^) [cf. 
Eq. (p^)] by the perturbative expression E°, Eq. ([ij), the difference contributing at sub-subleading order to the gap 
equation. In order to solve the gap equation, let us revert the analytic continuation to real energies in Eq. (^^, i.e., 
go is purely imaginary in the following. From Eqs. (|2|) and (|lj), the effect of including E° is to replace 



90 



90 



Ziqo) 



(15) 



in the quark propagator, where 



Ziqo)= 



(16) 



is the quark wave-function renormalization factor. Since we only want to consider the real part of the quark self- 
energy, we shall ignore the cut of the logarithm in Eq. (16) when performing the Matsubara sum in Eq. (O) by contour 



integration. In other words, we assume that the quark propagator has only simple poles in the complex go-plane, 
corresponding to the excitation energies of quasi-particles with infinite lifetime. This approximation is valid up to 
subleading order in the gap equation, because, as explained above, effects from a finite quasi-particle lifetime enter 
only at sub-subleading order. 
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A wave function renormalization of the form (16) is known from non-relativistic systems JTst , where it leads to 
non-Fermi liquid behavior. In relativistic systems, non-Fermi liquid behavior has been recently studied in great detail 
by Boyanovsky and de Vega [ p^ . 

After these introductory remarks, we may immediately proceed to Eq. (3.3) of ||ll| or Eq. (32) of This equation 
determines the spin-zero gap in a two-flavor color superconductor in pure Coulomb gauge. With the replacement ( p^ ) 
it reads 



*(-»■) 



+AtiK-Q) 
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MK~Q) 



1+kq 



k-q l + k-q(fc-g)2 



(17) 



where we neglected the contribution of anti-particles. The next step is to perform the Matsubara sum over go- We 
use spectral representations for the propagators, as in jsj. The only difference to the calculation of j|] is that the 
poles of the fermion propagator are shifted. To leading order, they are now given by 



go ^ ±^(eg) Eg = ±e. 



(18) 



The rest of the calculation is straightforward. We also take the external quark energy fco to be on the new quasi- 
particle mass-shell, fco = ■^(efc)efc = Cfc- Then, in analogy to Eq. (3.4) of Q and Eq. (72) of ||[, the final result for 
the gap equation, including the quark self-energy, reads 



djq - A^) ^2f~ 

^ l£<3 





tanh ( 

2T 




(19) 



where b = 256Tr^\2 / {Nfg^)f^^ . Note that we have replaced the symbol b in Eq. (72) of ||] by b, because the definition 
of b, cf. Eqs. (H^), includes 6o, the value of which has yet to be determined. In Ref. [H, this distinction was not 
necessary, because there &o = 1. We abbreviated 4>k = (/!)(efc,k); 0g is defined similarly. 



III. SOLVING THE GAP EQUATION 

Let us now solve the gap equation ( pj| ) at zero temperature. In this case, the factor tanh [eq/(2T)] = 1. Moreover, 
to leading order we can make the replacements — ^ and ik ^ ffc in the logarithm In (i^lJ-^ /{Sg — ■ For similar 
reasons, Z{eq) ~ Z{eq). Following Ref. 0, we approximate 

and then introduce the variables [|| 

y-^mf^VV (21b) 

.* = -gln(^^y (21c) 

xo = 5 In (^^^ . (21d) 

Note that in contrast to we choose to include a factor g in the definition of these variables. Consequently, since 
(j)o Mexp(— 1/g), X* ^ 0(1) and xo ^ 0(g). Furthermore, x and y are of order 0(1) near and of order 0(g) away 
from the Fermi surface. 
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In analogy to Eqs. (84) and (85) of H, the gap equation and its derivatives read in these new variables 



{x)~x I dy{l - 2gy)(t){y) + dy y{l - 2gy)(t){y) , 



d(l)(x) 

dx 
d'^(j){x) 
dx"^ 



dy (1 - 2gy) <j>{y) 
1 - 2gx) <t){x) . 



(22a) 
(22b) 
(22c) 



In these equations, we neglected contributions of order 0{g^), for instance, a term g^hi{bfj,/M) in the wave function 



renormalization factor Z(eo 



In order to solve Eq. (22c), we replace x with a new variable z, 

z = -(2g)-2/3(i_2gx) , 
and obtain Airy's differential equation ]l7t , 

d^(biz) , , , 
^ = -Hz) ■ 

The solution 0(z) of Eq. ( p4| ) is a linear combination of the Airy functions Ai(z) and Bi(z) 

(j){z) = Ci Ai(z) + Ca Bi(z) . 



(23) 



(24) 



(25) 



In weak coupling, z is always negative, and the Airy functions and their first derivatives can be expressed in terms of 
modulus and phase, defined as 

Ai(z) = MdzDcos^dzl) , Bi(z) = MdzDsin^dzl) , 

■Bi(z) 



Af(|z|) = JAi'^(z)+Bi'^(z) , 9(\z\) = arctan , , , 

Ai'(z) = 7V(|z|)cos<^(|z|) , Bi'(z) = A^(|z|)sin(p(|z|) , 

"Bi'(z) 



N{\z\) = y'Ai'2(z) +Bi'2(z) , <^(|z|) = arctan 



Ai'(z) 



(26) 



At the Fermi surface, the value of the zero-temperature gap function is (j){z*) — (jjo and its derivative vanishes, 
d(j){z*)/dz = 0. Consequently, we obtain for the gap function 



M{\z\) sin [ip{\z*\)-0i\z\)] 
M{\z*\) sin[^(|z*|)-0(|z*|)] 



(27) 



In order to determine (f>Oj we use Eq. (22a) at the Fermi surface, z = z*, and substitute the integration variable y 
by u = -(2g)~2/3(i _ 2gy) to obtain 



z*)^ du[u + {2g)-^/^]u(j){u) , 



(28) 



where zq = —(2g) ^/^[l — 2gxo]. According to Eq. (|2^), we can replace u(f){u) with d^ (j)(u) / dv?' . Integrating by parts, 
this leads to the condition 

[zo + (25)-'/^] 4>'{z^) = 0(zo) . (29) 

Note that the above equation depends on z* through Eq. (|2^). It seems that Eq. (|2|) also depends on zq which is 
arbitrary and far from the Fermi surface. In weak coupling, however, the dependence on zq disappears, as we shall 
show in the following. We first rewrite the condition (E9[) as 



[zo + (2ff)-2/3] sin 



M{\zo\) 
N{\z,\) 



sm 



v{\z*\)-e{]z,\) 



(30) 
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In weak coupling, \z\ ^ (2g) ^/'^ ^ 1, and we may use the asymptotic formulas jlV 

\ hx-g \ + O(g^) 

0{\z\) 



24 



iV(|z|) 



1-1/2 



1 + 0(1^1-3) 



where we employed \z\ ~ {2g) ^/^(l — 25^), cf. Eq. (p^). We now expand Eq. ( ^ ) to order 0{g) and obtain 

2 



X = arctan 



g{l + a;*2) 



(31) 



(32) 



In weak coupling, the argument of the arctan is large, and we can expand the right-hand side to order 0{g) around 
7r/2. The result is 



TT 1 + X 



*2 



(33) 

To order 0(5), we can approximate x*^ ~ tt'^ /A on the right-hand side of Eq. (|3^), and using the definition of x* , Eq. 
(21c), we obtain the zero-temperature gap value at the Fermi surface 



ho = 2 b bQ ^ exp 



TT 

2i 



(34) 



where &q is given by Eq. (^, with c[ of Eq. (||). In conclusion, the effect of including the quark self-energy in the gap 
equation changes the value of &q from one, as in MM, to the value c'l given in Eq. m). 



IV. DETERMINING THE TRANSITION TEMPERATURE 



Within the gap equation approach, we can also determine the temperature Tc for the transition between the normal 
and the superconducting phase. We follow Ref. Q and consider the gap equation (22a) at the Fermi surface, restoring 
the factor tanh[e(t/)/(2T)] present at non-zero temperature. As in B, we assume that to leading order the shape of 
the gap function at non-zero temperature is the same as at zero temperature, and that only the overall magnitude 
changes with temperature, 4>{x,T) = (j){T)(j}{x ,Q) / (pQ. Then the gap equation reads 



dy y(l - 2gy) tanh 



2T (bo 



(35) 



We now separate the range of integration into two pieces, [0:0,2;*] — > [xo,^;^] + [xk,^*], where x^. = x* — .gln(2K) = 
gln[6/i/(K(/)o)]. The main contribution to the integral in Eq. ( |35| ) comes from the region of momenta away from the 
Fermi surface, [xqjXk]. In this region, e{y) ^ <f>o ^T, such that we can approximate the factor tanh[e(?;)/(2T)] ~ 1. 
By making use of Eq. (pj) and Eq. (p9|), the integral over the region [xqjXk] is evaluated as 



dy 2/(1 - 2gy) 



0(y,o) 



1 



Z„ + (2.9) -2/3 <j,\z^) 



= l--g\u2K + 0{f) 



(36) 



where = z* — {2gY/^g\Ti2K,. The last line is obtained by expanding the right-hand side of the second equality to 
order 0{g) around z* . Equation (^sl) becomes 
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r d2/2/(l-2OT)tanh^^M^ Jgln2^. + 0(52). (37) 

The integral on the left-hand side may now be computed to order 0{g). As in this amounts to approximating 
y :^ X* and 4>{y, O)/0o — 1- Furthermore, the correction from the quark self-energy can be neglected, 1 — 2gy ~ 1. In 
this way, we obtain Eq. (104) of consequently the BCS result (|^) remains valid to leading order in g, even when 
the quark self-energy is taken into account in the gap equation. With Eqs. (^) and (|3^), we thus conclude that our 
result for Tc is the same as that obtained in Ref. 111]. 



V. CONCLUSIONS 



In this paper, we have computed the spin-zero gap in a two-flavor color superconductor at zero temperature from a 
mean-field gap equation. In contrast to ear her studies @,|,|,Hll), 

we have included subleading contributions from the 
real part of the quark self-energy. We found that these contributions reduce the gap parameter at the Fermi surface 
by a factor 6q = exp[— (tt^ -|- 4)/8] ~ 0.177. We then computed the transition temperature Tc between the normal 
and superconducting phase and found that the BCS relation Tc ~ 0.57(/!)o remains valid to leading order i n g after 
including the corrections from the quark self-energy. Therefore, we obtain the same value for Tc as in Ref. [[l2| . 
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